Lasercooled radium monofluoride: A molecular all-in-one probe for new physics 
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The particular advantages of using diatomic molecule radium monofluoride (RaF) as a versatile 
molecular probe for physics beyond the Standard Model are highlighted, i) RaF was previously 
suggested as being potentially amenable to direct cooling with lasers. As shown in the present 
work, RaF's energetically lowest electronically excited state is of ^ff symmetry (in contrast to BaF), 
such that no low-lying state prevents efficient optical cooling cycles, ii) The effective electric 
field acting on the unpaired electron in the electronic ground state of RaF is estimated larger than in 
YbF, from which the best restrictions on the electron electric dipole moment (eEDM) were obtained 
experimentally, iii) Favourable crossings of spin-rotational levels of opposite parity in external 
magnetic fields exist, which are important for the measurement of the nuclear anapole moment of 
nuclei with a valence neutron. Thus, RaF appears currently as one of the most attractive candidates 
for investigation of parity-odd as well as simultaneously parity- and time-reversal-odd interactions 
in the realms of molecular physics. 



INTRODUCTION 

The principle advantage of using heavy-atom, polar, 
diatomic molecules for searching for space parity vi- 
olation interactions (P-odd interactions) and simulta- 
neously space parity and time reversal violating forces 
{V, T-odd forces) is known for more than 30 years. Nev- 
ertheless, only the latest generation of molecular experi- 
ments finally surpassed their atomic competitors in sen- 
sitivity to one of the most important V, T-odd proper- 
ties of elementary particles, namely the permanent elec- 
tric dipole moment of an electron (eEDM) Whereas 
molecules arc ideally tailored to create favourable, well- 
defined fields at the heavy nucleus, the complexity of 
measurements with molecules is typically connected with 
various systematic effects which can mimic V-odd corre- 
lations. This calls for an active identification of promising 
molecular candidates that allow different approaches to 
suppress systematic effects. In the search for an eEDM, 
experimentally oriented research groups are currently fo- 
cussing on YbF 1], PbO 0, ThO^, WC 0, PbF @ 
molecules and HfH+ molecular ion [6| . In all these exper- 
iments, high-quality electronic structure calculations are 
crucial both for the preparation stage of the experiment 
and for subsequent interpretation of experimental data 
obtained 

Another set of pivotal molecular experiments is con- 
nected with attempts to measure the nuclear anapole 
moment, a V-odd electromagnetic form- factor appearing 
in / > nuclei due to V-odd nuclear forces. The only 
nucleus, in which the anapole moment was successfully 
determined, is ^^^Cs in experiment with a vapour of Cs 
atoms . The result of this experiment is apparently in 
disagreement with the earlier measurement on Tl atom 
[^, ■ Currently molecular experiments are under de- 
velopment in Yale on BaF ^ and in Groningen on SrF 

Recently, we identified the open-shell diatomic 
molecule RaF as an exceptionally suitable candidate for 



nuclear anapole moment measurements, having on the 
one hand a high enhancement factor for nuclear spin- 
dependent weak interaction and on the other hand offer- 
ing potential for direct cooling with lasers. In the present 
work we demonstrate that RaF presents unique possibili- 
ties for measurements of T'-odd and V, T-odd effects due 
to favourable combinations of peculiarities of molecular 
electronic structure and nuclear structure of isotopes of 
radium. 



DIRECT COOLING OF MOLECULES WITH 
LASERS 



We identified earlier a set of requirements on molec- 
ular electronic structure for molecules being suitable for 
direct cooling with lasers [lH. Monofluorides of group 
II elements (e.g. BaF and RaF) belong to the first class 
of molecules with highly diagonal Frank-Condon matrix. 
One problem emerging, however, in lasercooling of BaF 
molecule (isovalent to RaF) is connected with the exis- 
tence of a metastable level, lying energetically below 
the level involved in the optical cooling loop. Our 
previous electron correlation calculations of the spectro- 
scopic parameters in RaF indicated that the energeti- 
cally lowest electronically excited level is ^11. We took 
now larger atomic basis sets and active spaces of vir- 
tual molecular orbitals to investigate the stability of the 
ordering of electronic levels. The results are given in Ta- 
blelH These show that even considerable alterations in the 
parameters of the Fock-space relativistic coupled cluster 
(FS-RCC as implemented in Dirac [3]) calculations do 
not change the ordering of levels in RaF. On the other 
hand FS-RCC, calculations of BaF with the basis set of 
similar quality as for RaF and as large active spaces (see 
supplementary material) confirm the first excited elec- 
tronic level in BaF to be ^A. Based on comparison to 
experimental BaF results we estimate the accuracy of 
To calculation for RaF to be within 1200 cm~^ (without 
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changing the ordering of the levels), Rc within about 0.1 
ao and ojc about 60 cm"'^. 



LEVEL CROSSING IN MAGNETIC FIELD AND 
SENSITIVITY TO THE NUCLEAR ANAPOLE 
MOMENT 

■p-odd effects in diatomic molecules can be greatly 
enhanced by shifting levels of opposite parity to near- 
crossing with the help of external magnetic fields 17 [. 
This idea is exploited in for measurement attempts 



of the nuclear anapole moment in BaF. One of the main 
problems in the suggested approach is to create highly ho- 
mogeneous magnetic fields in large volumes. Favourable 
values of magnetic fields required to tune spin-rotational 
levels to near crossing would be below 10 kG (1 T), as 
creation of larger magnetic flux densities require special 
effort, for instance superconducting magnets. To esti- 
mate if fields with \B\ < IT suffice to create near level 
crossings in RaF we calculated Zeeman splittings for spin- 
rotational levels (See Fig. I and Fig. II in supplementary 
materials). Matrix elements of the spin-rotational Hamil- 
tonian in magnetic field were implemented as in Ref . [l8| . 
The following parameters of the spin-rotational Hamilto- 
nian were employed: rotational constant = 5689 MHz 
(as calculated from equilibrium structure), ratio of spin- 
doubling constant A to 2B, A/2B = 0.97 (as calcu- 
lated within a four-component Dirac-Hartree-Fock ap- 
proach), components of the hyperfine tensor for ^^'^Ra 
nucleus A\\ = -15100 MHz and Aj_ = -14800 MHz (cal- 
culated with the ZORA approach, values were already 
used for scaling in Ref. [l9|, but not explicitly reported 
therein) and components of G-tensor Gy = 1.993 and 
G_L = 1.961 (crude estimate based on results for HgH). 
According to our calculation, the first crossing of levels 
of opposite parity takes place at about 3 kG for levels 
with the projection Fz of the total angular momentum 
on the direction of the magnetic field being —3/2. A few 
more crossings take place in fields up to 10 kG for levels 
with Fz = ^1/2, thus providing additional freedom for 
choosing the optimal experimental parameters. 

To estimate the lowest possible flux of RaF molecules, 
which allows to measure the anapole moment in RaF, 
one needs to find the ratio with W being the experi- 
mentally measured signal proportional to matrix element 
of the nuclear spin-dependent weak interaction W^, and 
AW its experimental uncertainty. Of course the condi- 
tion for meaningful measurement is < 1 . We assume 
that the experimental scheme employed for anapole mo- 
ment measurement is analogous to the one suggested in 
According to Q the maximal value of 4^ is 
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AW 



W 2^/2NotW 



(1) 



where F is detected molecular flux and r is the to- 
tal measurement time) and t the interaction time be- 
tween molecule and external fields. Thus one gets F > 
l/(8VF^t^r). The time for molecular trapping can reach a 
few seconds [2^, so let i = 1 s, t = 1 h = 3600 s. To esti- 
mate W we just scale the W value for Ba given in Ref. , 

as VF^a/p^Ra ^ ^Ba/^Ra^ ^.^^j^g |^Ra ^ IQW^"" = 

50 Hz. The least required fiux of RaF is roughly then 
F 1/(8 • 2500 • 1 • 3600) = 1.4 • 10"^ s^^. In practice 
one might expect trapping and detection of at least one 
molecule during an experiment time of t = 1 h (see be- 
low), which corresponds to the flux 2.8 • 10^'* s^^, which 
is a few orders of magnitude higher than the least re- 
quired fiux of RaF. In practice this implies that it should 
be possible to perform successful measurements with sig- 
nals from single trapped RaF molecules. 



EFFECTIVE ELECTRIC FIELD ACTING ON 
THE UNPAIRED ELECTRON IN RAF 

One of the most important parameters in molecular 
experiments on eEDM is the effective electric field acting 
on the unpaired electron in the electronic ground state 
state of the molecule of interest. This field, however, 
can not be measured directly in experiment, but is pre- 
dicted from quantum chemical calculations. We estimate 
here the effective electric field acting on the unpaired 
electron in RaF by using relations between matrix ele- 
ments of different V-odd and VjT-odd operators as it 
has been done in Ref. ^21] by Kozlov and extended re- 
cently in R ef. |22t . According to the semiempirical model 
of Kozlov [21] the relation between the parameter Ws 
of the 7^,T-odd term and the parameter W^ of the V- 
odd term in the effective spin-rotational Hamiltonian is 
Ws/Wa — Z3'-f/{2'-f + 1), where Z is the nuclear charge 
number and 7 — y^l — (aZ)^. 

On the other hand, matrix elements Mspt of the scalar 
7',T-odd interaction and matrix elements Medm of the 
coupling between eEDM and inner molecular electric field 
can be expressed as [22]: 



Mspt/Gsp 

MEDu/dc 



'^^^-aAsApPo{A)Rj,{A){l - 0.56Z^a^()2) 



3^2 



7(472-1)1/, 



372 372" 



(3) 



In the equation above Gp is Fermi's constant, A the 
atomic mass number, i?N(A) = l.2A^^^ fm the effective 
nuclear radius and po{A) — An/ {3Rtss{A)^). The parame- 
ters As and Ap are normalisation factors of the electronic 
wavefunction and can be expressed in the semiclassical 
approximation as 22l |: 

(l+7)(2Zi?N(A))^ 



As^ 



with A^o being the total number of molecules (A'o = Ft, 



(1 - 0.375Z2a2)y^r(27 -I- 1 
Z{2ZR^{A))-' 



3/2 



(1 - 0.375Z2a2)r(27- 



. N 3/2 ' 



(4) 



(5) 



3 



TABLE I. Estimated molecular spectroscopic parameters for the electronic ground state and the first electronically 

excited state from FS-CCSD calculations of ^^^RaF and ^''^BaF. 17 electrons are correlated in both cases: 2pi/2,3/2 and 2si/2 
on fluorine, 5si/2, 5pi/2,3/2 and 6si/2 and Ba and 6si/2, 6pi/2,3/2 and 7si/2 an Ra. The results are given for two different 
4-component all-electron uncontracted basis sets on Ra and Ba atom: Dyall's relativistic basis set and the RCC-ANO basis set 
by Roos et al. The active space in FS-CCSD was restricted by energy for Dyall's basis set up to 10 Hartree and for RCC-ANO 
up to 1000 Hartree (see supplementary material for computational details and references). 

Rc/ao LOc/cm^^ TL/cm"^ Dc/cm~^ 
RaF 
Dyall basis set 







4.24" 


428" 




32104" 


'ni/2 




4.24" 


432" 


14012" 


31302" 


'n3/2 


A3/2 


4.25 


410 


15967 




'A3/2 


n3/2 


4.25 


432 


16435 




'A5/2 




4.27 


419 


17086 




^Si/2 




4.26 


416 


18118 








RCC-ANO basis set 




^Si/2 




4.29 


431 




52425 


'ni/2 




4.29 


428 


13263 




'n3/2 


A3/2 


4.31 


415 


15020 




'A3/2 


113/2 


4.28 


431 


15405 




'A5/2 




4.30 


423 


15840 




^Si/2 




4.32 


419 


16682 










BaF 










RCC-ANO basis set 




^Si/2 




4.15 


456 




45006 


^A3/2 




4.21 


455 


10871 




'A5/2 




4.21 


455 


11304 




'ni/2 




4.19 


446 


11701 




^113/2 




4.19 


444 


12354 




^Si/2 




4.23 


455 


14220 





Experimental data 
(from Ref. [14] if not indicated otherwise) 



Sl/2 



4.09" 

4.13^= 
4.17^= 



469'' 
437 
437 
424 



46780 ± 726" 



10940 
11727 
13828 



a) Ref. [il; b) Ref. [H; c) Ref. [le 



with r(x) being the gamma- function and v the effective 
principal quantum number. Upon relating these quan- 
tities, one has to take into account different definitions 
of the parameters of the effective Hamiltonian and P, T- 
odd parameters {Ws and Wc) in '22^ and [21|] (see ^ for 
details). 

It is interesting to note that in experiments with one 
kind of molecules the scalar T', T-odd interaction is indis- 
tinguishable from the eEDM effect - they both contribute 
to V, T-odd electron paramagnetic resonance signal. One 
can, however, separate the contributions by taking data 
from experiments with different molecules (or molecule 
and atom) like it is proposed in ^2^] for eEDM in Tl and 
YbF. 



Taking into account the above relations one can easily 
estimate the effective electric field acting on the electron 
in the ground state of RaF. The results are given in 
Table |lTl The accuracy of such an estimate is not high, 
but sufficient to identify RaF as promising candidate for 
eEDM measurements. 

Another attractive feature of measurements with Ra 
nuclei is that there exist also a nuclear mechanism en- 
hancing T'-odd and V, T-odd effects in certain Ra iso- 
topes. According to [IBj and [2^ the Schiff moment in 
nuclei possessing octapole deformation is enhanced by 
about 10 to 100 times. The mechanism is similar to the 
one in diatomic and chiral molecules: enhancement is 
reached due to closeness of rotational levels of opposite 
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TABLE II. Estimates of the \ Ws \ parameter of P, T-odd scalar 
interaction (in kHz) and effective field on the electron |i5efi| 
(in units 10^"* Hz e"^ cm ) for open-shell diatomic molecules. 
For comparison the data of ab initio calculations of l-Ecff | arc 
provided. 

Z IVy.l/Hz") |W^,|/kHz \E,b\ \E,a\ 

(ab initio) 



BaF 


56 


1.9x10^ 




12 


1.8 


1.9') 


YbF 


70 


6.1x 10^ 




38 


4.1 


6.0") 


RaF 


88 


2.1x 10^ 


13x 


10^ 


9.5 




HgH 


80 


2. Ox 10^ 


19x 


10^ 


16 




CnH 


112 


3.1x 10* 


35 X 


10^ 


300 





a) Results of Ref. [Ig 

b) Results of Ref. [H 

c) Results of Ref. '2^ 



parity. As a result the estimated SchifF moments in ^^^Ra 
and ^^^Ra isotopes are equal to 300 and 400 (in units 
7y 10^ e fm'^, where rj is effective nucleon-nucleon V,T- 
odd force constant), respectively, whereas for ^^^Hg for 
example it is only —1.4 (data taken from Ref. [26|). 



PRODUCTION OF RAF 
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Besides the possible routes discussed in Ref. 
propose in Ref. [Ill to produce neutral RaF via RaF+, 
which is subsequently neutralised by charge exchange in 
collision with a suitably chosen collision gas or by in- 
teraction with surfaces that provide the adequate work 
function for an iso-enthalpic electron transfer. RaF"*" can 
in turn be formed in reactive collisions of radium ions 
with a suitable fluorine containing compound. 



CONCLUSION 

We demonstrated various special properties of RaF 
that render it a versatile molecular laboratory for study- 
ing wide range of physical phenomena, from laser cool- 
ing to the physics beyond the Standard Model. Unique 
combination of rovibronic, electronic and nuclear struc- 
ture features render RaF particularly attractive for fur- 
ther experimental study. The authors are grateful to D. 
DeMille, V. Flambaum, M. Kozlov and S. Hoekstra for 
discussion. 
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